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^ (57) Abstract: A fast, efficient and accurate pseudo 2-D inversion scheme (Fig. 4) for resistivity determination of an anisoUopic 
^ fomation (Fig. 3) uses dau from a tool (3DEX) (10) comprising three uansmitiers (101, 103, 105) and three corresponding re- 
Q ceivers ( 107, 109, 1 II) sampling the formation in a plurality of spatial directions. An initial model (311) of the fomiation including 
invasion zones is obtained using a conventional multifrequency and/or multi spacing logging tool. A pseudo 2-D inversion scheme 
Q combines an accurate full 2-D forward solution (313) of the synthetic responses of the earth model with a 1-D approximation (323) 
^ of the sensitivity matrix (317) of the horizontally layered anisotropic background model. The timesaving compared to a regular 2-D 
^ inveision scheme can be uemendous. 
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AN EFFICIENT AND ACCURATE PSEUDO 2-D INVERSION SCHEME FOR 



MULTICOMPONENT INDUCTION LOG DATA 

Berthold KriegshSuser, Steve Mc Williams, Otto Fanini; Liming Yu 

BACKGROUND OF THE INVENTION 
L Field of the InventioQ 

[0001] The invention is related generally to the field of interpretation of 
measurements made by well logging instruments for the purpose of determining the 
properties of earth formations. More specifically, the invention is related to a method 
for inversion of measurements made by multi-component induction or propagation 
sensors. 

2. Background of the Art 

[0002] Electromagnetic induction and wave propagation logging tools are commonly 
used for determination of electrical properties of formations surrounding a borehole. 
These logging tools give measurements of apparent resistivity (or conductivity) of the 
formation that when properly interpreted are diagnostic of the petrophysical 
properties of the formation and the fluids therein. 

[0003] Tlie physical principles of electromagnetic induction resistivity well logging 
are described, for example, in, H. G. Doll, Introduction to Induction Logging and 
Application to Logging of Wells Drilled with Oil Based Mud, Journal of Petroleum 
Technology, vol. 1, p. 148, Society of Petroleum Engineers, Richardson Tex. (1949). 
Many improvements and modifications to electromagnetic induction resistivity 
instruments have been devised since publication of the Doll reference, supra. 
Examples of such modifications and improvements can be found, for example, in U.S. 
Pat. No. 4,837,51 7; U.S. Pat. No. 5,157,605 issued to Chandler et al, and U.S. Pat. 
No. 5,452,761 issued to Beard cl al. 

[0004] A limitation to the electromagnetic induction resistivity well logging 
instruments known in the art is that they typically include transmitter coils and 
receiver coils wound such that the magnetic moments of these coils are substantially 



1 
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parallel only to the axis of the instrument. Eddy currents are induced in the earth 
formations from the magnetic field generated by the transmitter coil, and in the 
induction instruments known in the art these eddy currents tend to flow in ground 
loops which are substantially perpendicular to the axis of the instrument. Voltages are 
then induced in the receiver coils related to the magnitude of the eddy currents. 
Certain earth formations, however, consist of thin layers of electrically conductive 
materials interleaved with thin layers of substantially non-conductive material. The 
response of the typical electromagnetic induction resistivity well logging instrument 
will be largely dependent on the conductivity of the conductive layers when the layers 
are substantially parallel to the flow path of the eddy currents. The substantially non- 
conductive layers will contribute only a small amount to the overall response of the 
instrument and therefore their presence will typically be masked by the presence of 
the conductive layers. The non-conductive layers, however, are the ones which are 
typically hydrocarbon-bearing and are of the most interest to the instrument user. 
Some earth formations which might be of commercial interest therefore may be 
overlooked by interpreting a well log made using the electromagnetic induction 
resistivity well logging instruments known in the art. 

[0005] Baker Atlas and Shell International E&P jointly developed a new 
multicomponent induction logging tool, 3DEX to measure the electrical anisotropy of 
these sequences. This logging tool and its use is described in United States Patent 
6,147,496 to Strack et al. . The instrument comprises three mutually orthogonal 
transmitter-receiver configurations that provide all necessary data to compute 
horizontal and vertical resistivities of the formation. These resistivities may then be 
used in an integrated petrophysical analysis to provide an improved estimate of the 
laminar sand resistivity and corresponding net oil-in-place. 

(00081 United States Patent 5,999,883 issued to Gupta et al , (the "Gupta patent"), the 
contents of which are fully incorporated here by reference, discloses a method for 
determination of the horizontal and vertical conductivity of anisotropic earth 
formations. Electromagnetic induction signals induced by induction transmitters 
oriented along three mutually orthogonal axes are measured. One of the mutually 



wo 02/071099 PCT/US02/05797 

orthogonal axes is substantially parallel to a logging instrument axis. The 
electromagnetic induction signals are measured using first receivers each having a 
magnetic moment parallel to one of the orthogonal axes and using second receivers 
each having a magnetic moment perpendicular to a one of the orthogonal axes which 

S is also perpendicular to the instrument axis. A relative angle of rotation of the 
perpendicular one of the orthogonal axes is calculated from the receiver signals 
measured perpendicular to the instrument axis. An intermediate measurement tensor 
is calculated by rotating magnitudes of the receiver signals through a negative of the 
angle of rotation. A relative angle of inclination of one of the orthogonal axes, which 

10 is parallel to the axis of the instrument is calculated, from the rotated magnitudes, 
with respect to a direction of the vertical conductivity. The rotated magnitudes are 
rotated through a negative of the angle of inclination. Horizontal conductivity is 
calculated from the magnitudes of the receiver signals after the second step of 
rotation. An anisotropy parameter is calculated from the receiver signal magnitudes 

1 5 after the second step of rotation. Vertical conductivity is calculated from the 
horizontal conductivity and the anisotropy parameter. 

[0009] However, the new horizontal magnetic field responses, H„ and H„, that are 
sensitive to the vertical resistivity of the formation can suffer from strong borehole 
20 and near-zone effects. These effects increase with borehole size, borehole fluid 
conductivity, and invasion depth. Co-pending United States Application Ser. No. 
09/676,097 by Kriegshduser et al. describes a method for applying shoulder bed 
corrections to this data. However, in large borehole sizes and conductive mud 
systems, these corrections cannot completely eliminate the near-zone effects. 

25 

lOOlOJ However, a rigorous 2-D inversion is too time-intensive and prohibitive for 
typical logging applications, because the 2-D forward model response is time 
consuming and also because inversion schemes typically require the determination of 
a Jacobian matrix defining the sensitivity of each of the measurements to every one of 
30 the parameters in the model. 

(0011 J Tabarovsky and Rabinovich (US Patent 5703773) teach a computationally fast 

3 
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method for 2-D inversion of induction logging data. The method includes skin effect 
correcting the responses of the receivers by extrapolating the receiver responses to 
zero frequency. A model is generated of the media surrounding said instrument. 
Conductivities of elements in the model are then adjusted so that a measure of misfit 
S between the skin-effect corrected receiver responses and simulated receiver responses 
based on the model is minimized. The geometry of the model is then adjusted so that 
the measure of misfit between the skin-effect corrected receiver responses and the 
simulated receiver responses based on the model is further minimized. In a preferred 
embodiment of the invention, the step of adjusting the geometry includes minimizing 
10 the measure of misfit between the simulated responses and the receiver responses 

from selected ones of the receivers closely spaced to the transmitter. Numbers of and 
positions of radial boundaries are then determined by minimizing the measure of 
misfit for all the receiver responses. However, Tabarovsky does not address the 
problem of inversion of multi-component data. 

15 

[0012] There is a need for a method of inversion of multicomponent induction 
logging data that gives reasonably accurate results without using an inordinate 
amount of computer time. The present invention satisfies this need. 

20 SUMMARY OF THE INVENTION 

[0013] An electromagnetic logging tool having a plurality of transmitters and 
receivers (3DEX) is used to obtain multicomponent measurements indicative of 
horizontal and vertical resistivities of subsurface formations. A model for the 

25 horizontal resistivity, length of the invasion zone and resistivity of the invasion zone 
may be obtained from High Definition Induction Logging (HDIL) tools. Such an 
induction logging tool comprises transmitter and receiver coils with axes parallel to 
the tool axis: measurements are made at multiple frequencies and/or with multiple 
transmitter-receiver spacings. An example of such a tool is given in U.S. Patent 

30 5,452,761 to Beard et al. An initial model is defined that includes the obtained model 
and vertical resistivities for the formations. In one embodiment of the invention, the 
vertical resistivities for the initial model are set equal to the horizontal resistivities. In 

4 
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FIG. 6 shows a comparison of inversion results using a 1-D inversion and the P2D 
inversion method of the present invention with a model. 

DETAILED DESCRIPTION OF THE INVENTION 

[0017] Referring now to Fig. 1, an electromagnetic induction well logging instrument 
10 is shown disposed in a wcUbore 2 drilled through earth formations. The earth 
formations are shown generally at 4. The instrument 10 can be lowered into and 
withdrawn from the wellbore 2 by means of an armored electrical cable 6 or similar 
conveyance known in the art. The instrument 10 can be assembled from three 
subsections: an auxiliary electronics unit 14 disposed at one end of the instrument 10; 
a coil mandrel unit 8 attached to the auxiliary electronics unit 14; and a 
receiver/signal processing/telemetry electronics unit 12 attached to the other end of 
the coil mandrel unit 8, this unit 12 typically being attached to the cable 6. 

[0018] The coil mandrel unit 8 includes induction transmitter and receiver coils, as 
will be further explained, for inducing electromagnetic fields in the earth formations 4 
and for receiving voltage signals induced by eddy currents flowing in the earth 
formations 4 as a result of the electromagnetic fields induced therein. 

(0019] The auxiliary electronics unit 14 can include a signal generator and power 
amplifiers (not shown) to cause alternating currents of selected frequencies to flow 
through transmitter coils in the coil mandrel unit 8. 

[0020] The receiver/signal processing/telemetiy electronics unit 12 can include 
receiver circuits (not shown) for detecting voltages induced in receiver coils in the 
coil mandrel unit 8, and circuits for processing these received voltages (not shovm) 
into signals representative of the conductivities of various layers, shown as 4A 
through 4F of the earth formations 4. As a matter of convenience the receiver/signal 
processing/telemetry electronics unit 12 can include signal telemetry to transmit the 
conductivity- related signals to the earth's surface along the cable 6 for further 
processing, or alternatively can store the conductivity related signals in an appropriate 
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recording device (not shown) for processing after the instrument 10 is withdrawn 
from the wellbore 2. 

(0021) Referring to Fig. 2, the configuration of transmitter and receiver coils in a 
5 preferred embodiment of the 3DExpiorer ™ induction logging instrument of Baker 
Hughes is shown. Three orthogonal transmitters 101, 103 and 105 that are referred to 
as the Tx, T^, and Ty transmitters are shown (the z- axis is the longitudinal axis of the 
tool). Corresponding to the transmitters 101, 103 and 105 are associated receivers 
107, 109 and 1 11, referred to as the R^, R^, and Ry receivers, for measuring the 
1 0 corresponding magnetic fields. In a preferred mode of operation of the tool, the 

Hyy /L, and H^^ components are measured, though other components may also 
be used. 

[0022] Fig. 3 is a schematic illustration of the model used in the present invention. 

1 5 The subsurface of the earth is characteri2;ed by a plurality of layers 201a, 201b, . • 
.201 i. The layers have thicknesses denoted by /r,, hy . . A,^. The horizontal and 
vertical resistivities in the layers are denoted by . . .R^, and R^j, R,^, . . ,R^, 

respectively. Note that the top and bottom layers are semi-infinite in the model. The 
borehole is indicated by 202 and associated with each of the layers are invaded zones 

20 in the vicinity of the borehole wherein borehole fluid has invaded the formation and 
altered its properties so that the electrical properties are not the same as in the 
uninvaded portion of the formation. The invaded zones have lengths Z^,, Z^^^, . . ,L^. 
extending away from the borehole. The resistivities in the invaded zones are altered 
to values R^oj, R^^^ • • -^jor I" the embodiment of the invention discussed here, the 

25 invaded zones are assumed to be isotropic while an alternate embodiment of the 
invention includes invaded zones that are anisotropic, i.e., they have different 
horizontal and vertical resistivities. The assumption of an isotropic invasion zone is 
reasonable because in the case that the borehole fluid is conductive and invades a 
laminated sand/shale layer, then the pore fluid of the sand laminae is filled with 

30 conductive borehole mud fluid. Hence, the sand laminae become as conductive as the 
shale laminae, making this invaded zone isotropic. If the borehole fluid is resistive 
and invades the sand/shale layer, then the resistive pore fluid of the sand laminae is 

7 
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replaced by resistive borehole fluid. Hence, there is no significant resistivity contrast 
between the invasion zone and the anisotropic formation layer. 

[0023] The observed data D may be defined as a function of the model parameters m 
S as 

D^f{m) 

or, in matrix form» 







' Aim)' 


d. 




Aim) 






A.Am) 


.d, . 




. fii^) . 



is the vector of observations 



ID 



where 



W = [m, nty W5 m^^^ mf^^ m^^^ ^i^t ^i] = 

(1) 

15 is a model vector comprising layer thicknesses, horizontal and vertical resistivities of 
the layers, length of the invasion zones in each of the layers and resistivity of the 
invasion zone in each of the layers. There are k layers in all, so that for the preferred 
embodiment of the invention where the invaded zone is characterized by a length of 
the invaded zone and a single resistivity, w is a vector of length 5k'2 For a full 2-D 

20 inversion, the data vector could comprise all the components of measurements made 
with the SDEX^*^ tool. In a preferred embodiment of the present invention, the 
component and at least one of (i) the component, (ii) the component, and, (iii) 
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an average of the ^ and 
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compared w,* ao^ j.,, n,.,^^^ ^^T^ ^ «f model 313 is 
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the entire set of observations using the entire 2-D model vector as defined above with 
reference to Fig. 3. This is v^ould be done by iteratively updating m to minimizie 

0,=||D^*'«/(m)f (2) 

5 

where D^' are the observations and Oj is a data objective function. 

(0027J In Fig. 4, the full 2-D inversion would involve the minimization of eq. (2), 
using the full model vector m from 321, and determining the full sensitivity matrix of 
10 to all the parameters in the model. This requires determination and inversion of a 

Jacobian matrix of sensitivities 



15 



dd 



d m J 



, which is a matrix comprising / rows and 5k -2 columns. In the 



present invention, a partial Jacobian matrix of sensitivities is used vnth respect to a 
subset m'of model parameters* 



[0028] The subset m 'of the model m comprises only the formation parameters, i.e., 
the layer thicknesses and resistivities, and the remaining components of the model m 
are assumed to be knovm. The invasion zone parameters, Rxo and Lxo, can for 
example, be derived from HDIL data, or any other a priori information. However, the 
20 inversion process does not alter these parameters of the invasion zone. In other 

words, the iterative procedure only updates the model m ' Note, however, that the 
forward model solution 313 is still a full 2-D model including the invasion zone 
effects. 



25 [0029] As a further approximation, the sensitivity matrix is determined for only a 
subset of the measurements 323 using only the and at least one of (i) the H„ 
component, (ii) the Hyy component, and, (iii) the average of the Hxx and Hyy 
components of the data. The sensitivity matrix, however, is determined for an 
objective function defined for a layered 1-D model of a horizontally layered 

10 
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formation. 



O' =Id<^' - Am')f (3) 

By using this approximation of the sensitivity matrix, computations are greatly 
S speeded up. 

EXAMPLE 

[0030J Next, an illustrative example of the use of the pseudo 2-D inversion scheme 
10 on a simple but yet realistic and typical 2-D model is shown. The formation 

parameters of the layered anisotropic background model are shown in Table 1 . The 
borehole diameter is 12.25 inches, and the borehole mud resistivity is O.I Ohm-m. 
The syntiietic example shown is typical for well logging applications, when seawater 
is used as drilling fluid. 

15 



TABLE I: RESISTIVITY MODEL PARAMETERS 



Layer boundaries 


R.(Q-m) 


(Q-m) 


<0 


0.8 


1.2 


36893 


11 
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36955 


1.5 


7.9 


37048 


10 


12 


37080 


1 


1.3 


37147 


15 


17 


14-16 


1 


1.3 


16-20 


12 


15 


20-22 


2.4 


10.7 


22-24 


1 


4.3 


>24 


0.9 


0.9 



10031) The parameter resolution can be studied qualitatively witii ♦back-transformed 
singular values' (BTSV) (Jupp and Vozoff, 1977; Hdrdt, 1992). To define tht BTSV, 
we begin with the singular value decomposition of the Jacobian matrix J, J=USV\ 
where U is the orthonormal matrix of data eigenvectors, V is tiie orthonormal 
parameter eigenvector matrix and S is the diagonal mau-ix containing Uie singular 
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values of J^J (e.g., Jupp and Vozoff, 1974). 

I0032J Now. we introduce the diagonal matrix K, composed of the normalized 
singular values (A; = Sj/s(l),j = 7.2. . . m) as diagonal elements. S(l) is the largest 
singular value; m is the number of parameters. The norm of each column of the 
matrix VK is called the back-transformed singular value (BTSV), and is considered to 
be a measure of the importance of the contribution of parameter to reduce the data 
misfit in the inverse process. BTSV varies between zero and one, where the value is 
directly proportional to the importance of the parameter. 

(0033] Fig. 5 shows the corresponding BTSV for both the 1-D and 2-D formation 
models. The bar chart clearly shows that the BTSV are almost identical for the 1-D 
and the 2-D models. This similarity indicates that a 1-D Jacobian can be used in the 
pseudo 2-D inversion process and that the formation parameter resolution should be 
almost identical compared to the full 2-D inversion. 

[0034) These synthetic responses are then inverted using a 1-D inversion scheme and 
the results compared with the final iterate of the pseudo 2-D inversion scheme (Fig. 
6). In the inverse process we only inverted for horizontal and vertical resistivities of 
the formation layers. The 1-D inversion scheme assumes no borehole and invasion 
zones in the synthetic formation model. Layer boundaries and borehole parameters 
are known and fixed during inversion. This assumption is usually valid in logging 
applications. Fig. 6 shows a comparison of the true formation resistivity parameters, 
the recovered model parameters using a 1-D inversion scheme, and the pseudo 2-D 
fmal iterates (P2D). 

(00351 The horizontal resistivity values for the 1-D inversion and the P2D inversion 
are substantially the same as the true value 401. The vertical resistivity obtained by 
the P2D is substantially identical to the true value 411 but the 1-D inversion 413 is 
clearly erroneous at certain depths, so that the determined resistivity ratio from the 
P2D is correct 421 while the determined resistivity ratio from the 1-D inversion 423 
deviates from the correct value. 
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(0036] For this simple model, the computation time for convergence of the P2D 
method on a SUN Ultra 5 workstation is comparable to that of 1-D inversion for the 
horizontal resistivity while the 2-D inversion took approximately 5 times as long. For 
5 the vertical resistivity, the P2D method took about 1 2 times as long as the 1 -D 

inversion while the 2-D inversion took about 7 times as long as the P2D inversion. 
There is thus a significant saving in computation as a result of using the method of the 
present invention over a full 2-D inversion with little loss in accuracy. 

1 0 (0037J The present invention has been discussed above with respect to measurements 
made by a transverse induction logging tool conveyed on a wireline. This is not 
intended to be a limitation and the method is equally applicable to measurements 
made using a comparable tool conveyed on a measurement-while-drilling (MWD) 
assembly or on coiled tubing. 

15 

(00381 While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to those skilled in the art. It is 
intended that all variations within the scope and spirit of the appended claims be 
embraced by the foregoing disclosure. 

20 , , 
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1 1 . A method of determining a parameter of interest of subsurface formations 

2 surrounding a borehole, the method comprising: 

3 (a) obtaining a plurality of measurements indicative of a horizontal and 

4 vertical resistivity of the formation; 

5 (b) obtaining a plurality of layer intervals for the subsurface; 

6 (c) defining an initial model for the subsurface; 

7 (d) determining expected values for said plurality of measurements based 

8 on said model; 

9 (e) determining a difference between said expected values and actual 
1 0 values of said plurality of said measurements; 

^ ' (0 determining a sensitivity matrix relating a subset of said plurality of 

1 2 measurements to a subset of said model; and 

13 (g) iteratively updating said model by using said subset of the plurality of 

1 4 measurements and said sensitivity matrix, and reducing said 

15 difference. 

1 2. The method of claim 1 wherein obtaining said plurality of measurements 

2 further comprises conveying an electromagnetic logging tool into the borehole 

3 and making said measurements using a plurality of combinations of 

4 transmitters and receivers on the logging tooL 

1 3. The method of claim I wherein said electromagnetic logging tool comprises a 

2 first transmitter having an axis substantially parallel to an axis of the borehole. 

3 and second and third transmitters having axes orthogonal to said axis of the 

4 borehole and to each other. 

1 4. The method of claim 1 wherein obtaining said plurality of layer intervals 

2 further comprises using a HDIL logging tool. 
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1 5. The method of claim 1 wherein defining said initial model further comprises 

2 using an HDIL logging tool. 

1 6. The method of claim 1 wherein determining expected values of said 

2 measurements further comprises using a 2-D modeling program. 

1 7. The method of claim 1 wherein said subset of measurements comprises at 

2 least one of (i) a //„ component, (ii) a //^ component, and, (iii) an average of 

3 aH„ component and a H^y component. 

1 8. The method of claim 1 wherein said subset of said model comprises horizontal 

2 and vertical resistivities associated with said plurality of layers. 

1 9. The method of claim 1 wherein iteratively updating said model further 

2 comprises using the Marquardt-Levenberg algorithm. 

1 1 0. The method of claim 1 wherein said model comprises, for each of said 

2 - r plurality of layers, (i) a thickness, (ii) a horizontal resistivity, (iii) a vertical 

3 resistivity, (iv) a length of an invaded zone, and, (v) a resistivity of an invaded 

4 zone. 

1 II. A method of determining a parameter of interest of subsurface formations 

2 surrounding a borehole, the method comprising: 

3 (a) conveying an electromagnetic logging tool into the borehole and using 

4 at least one electromagnetic transmitter and one electromagnetic 

5 receiver on the tool for obtaining a plurality of measurements 

6 indicative of a horizontal and vertical resistivity of the formations; 

7 (b) obtaining a plurality of layer intervals for the subsurface, each of said 

8 layer intervals having a thickness; 

9 (c) defining an initial model for the subsurface, said subsurface model 

1 0 comprising, for.each of said plurality of layer intervals, a horizontal 

1 1 resistivity, a vertical resistivity, a length of an invasion zone, and a 

15 
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resistivity for said invasion zone; 


13 




determining expected values for said plurality of measurements based 


14 




on said model; 


IS 


(e) 


determining a difference between said expected values and actual 


16 




values of said plurality of said measurements; 


17 


(0 . 


determining a sensitivity matrix relating a subset of said plurality of 


18 




measurements to a subset of said model; and 


19 


(8) 


iteratively updating said model by using said subset of the plurality of 


20 
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21 




difference. 


1 


12. A method of determining a parameter of interest of subsurface formations 


2 


surrounding a borehole, said formation including invaded zones, the method 


3 


comprising: 


4 


(a) 


obtaining a plurality of measurements indicative of a horizontal and 


5 




vertical resistivity of the formation; 


6 


(b) 


obtaining a plurality of layer intervals for the subsurface; 


7 


(c) 


defining an initial model for the subsurface; 


8 


- (d) 


determining expected values for said plurality of measurements based 


9 




on said model; 


10 


(e) 


determining a difference between said expected values and actual 


11 




values of said plurality of said measurements; 


12 


(0 


determining a sensitivity matrix relating a subset of said plurality of 


13 




measurements to a subset of said model wherein said sensitivity matrix 


14 




does not depend upon parameters of said invaded zones; and 


IS 


(g) 


iteratively updating said model by using said subset of the plurality of 


16 




measurements and said sensitivity matrix, and reducing said 


17 




difference. 



16 
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